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I)  ABSTRACT  - -  ' 

The  purpose  of  the  theoretical  effort  at  Princeton  is  to  obtain  as  complete  a 
theoretical  understanding  as  possible  of  not  only  the  important  physical  processes 
which  Ic'd  to  the  experimentally  observed  effects  in  the  current  ionospheric  modifi¬ 
cation  program  (being  undertaken  at  the  Arecibo  Observatory  and  at  Boulder,  Colorado) , 
but  also  the  scaling  laws  which  predict  the  results  of  the  future  modification  experi¬ 
ments  (e.g.,  increased  high-frequency  modifier  power  fluxes). 

Extensive  diagnostics  reveal  distinct  collective  phenomena  on  several  scale 
lengths.  Our  understanding  of  short-scale  turbulence  (wavelengths  of  10^-10^  cm) 
allows  us  to  predict  quantitatively  the  level  of  density  fluctuations  which  can  be 
produced  in  a  controlled  manner  in  the  ionosphere.  We  show  that  the  presence  of  in¬ 
tense  modifier  radiation  will  cause  plasma  (parametric)  instabilities  in  the  ionosphere 
which  result  in  short-scale  turbulence  involving  plasma  waves  and  ion  acoustic  waves, 
and  lead  to  an  increased  absorption  of  the  modifier  power.  The  density  fluctuations 
associated  with  the  turbulence  can  scatter  radar  waves  as  evidenced  by  tiie  enhanced 
signals  seen  by  the  Arecibo  430  MHz  diagnostics.  Our  theory  has  described  both  the 
intensity  and  frequency  spectrum  of  the  Arecibo  radar  scatter  diagnostics.  The 
theory  provides  «.  good  description  of  the  frequency  spectrum  of  the  short-scaie  density 
fluctuations,  especially  the  plasma  lin**  portion. 

Our  work  on  long-scale  turbulence  (wavelengths  greater  than  an  ion  Larmor  radius 
s  10  m)  concerns  a  theory  for  artificial  Spread-F,  based  on  a  model  of  a  radiation  - 
induced,  self-focussing  thermal  instability.  The  predicted  growth  times  are  in  agree¬ 
ment  with  the  observed  fomation  times  (on  the  order  of  1-10  seconds)  of  magnetic 
field  aligned  density  structures. 
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ABSTRACT 


,\ 

Collective  plasma  phenomena  induced  in  the  ionospheric  F-region 
by  intense  high-frequency  radiation  are  examined  on  several  time  and 
length  scales.  The  short  wavelengh,  high-frequency  turbulence  asso¬ 
ciated  with  the  observed  high-frequency  enhanced  electron  density 
fluctuations  at  the  electron  plasma  frequency  is  shown  to  result 
from  parametric  excitation  of  plasma  and  acoustic  waves  by  the  high- 
frequency  radiation.  The  wave  kinetic  equation  describing  the  non¬ 
linear  evolution  or  the  turbulence  is  presented.  The  time  asymptotic 
spectral  intensity  in  calculated  and  a  formula  for  the  total  energy 
absorption  from  the  incident  radiation  via  nonlinear  collective 
effects  is  obtained.  A  radiation  induced  self— focussing  instability 
is  found  on  a  length  scale  somewhat  larger  than  an  ion  Larmor  radius, 
and  is  offered  as  a  possible  explanation  for  the  occurrence  of  magne¬ 
tic  field  aligned  density  filaments,  which  are  associated  with  the 
phenomenon  of  "artificial  Spread-F."  The  calculated  growth  rates 
are  in  agreement  with  the  time  delay  for  appearance  of  spread-F  after 
modifier  turn  on. 
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I .  SUMMARY 


> 


The  purpose  of  the  theoretical  effort  at  Princeton  is  to  ob¬ 
tain  as  complete  a  theoretical  understanding  as  possible  of  not 
only  the  important  physical  processes  which  lead  to  the  experi¬ 
mentally  observed  effects  in  the  current  ionospheri  ;  modification 
program  (being  undertaken  at  the  Arecibo  Observatory  and  at 
Boulder,  Colorado) ,  but  also  the  scaling  laws  which  predict  the 
results  of  the  future  modification  experiments  (e.g.,  increased 
high-frequency  modifier  power  fluxes) . 

Extensive  diagnostics  reveal  distinct  collective  phenomena 
on  several  scale  lengths.  The  shortest  scale  phenomena  (101  - 
10^  cm)  are  enhanced  high-frequency  electron  plasma  oscillations 
concomittant  with  an  increase  in  the  level  of  low-frequency  (on 
the  order  of  the  ion  acoustic  frequency)  oscillations  at  the  same 
wavevector.  At  substantially  longer  wavelengths  (0.1-1  km)  density 
striations  with  wavevectors  nearly  orthogonal  to  the  geomagnetic 
field  are  observed. 

Our  understanding  of  short-scale  turbulence  allows  us  to  pre¬ 
dict  quantitatively  the  level  of  density  fluctuations  which  can  be 
produced  in  a  controlled  manner  in  the  ionosphere.  The  controlled 
production  of  short-scale  density  fluctuations  may  have  important 
technological  applications.  Specifically,  our  work  has  shown  that 
the  presence  of  intense  modifier  radiation  will  cause  plasma 
(parametric)  instabilities  in  the  ionosphere.  The  plasma  insta¬ 
bilities  result  in  short-scale  turbulence  involving  plasma  waves 
and  an  increased  absorption  of  the  modifier  power.  "Parametric 
instability"  means  that  electromagnetic  modifier  power  spontaneously 
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converts  itself  into  a  combination  of  electronic  plasma  waves 
and  acoustic  waves.  The  density  fluctuations  associated  with  the 
plasma  waves  and  acoustic  waves  can  scatter  radar  waves  as  evi¬ 
denced  by  the  enhanced  signals  seen  by  the  Arecibo  430  MHz  diag¬ 
nostics.  Our  theory  has  described  both  the  intensity  and  frequency 
spectrum  of  the  Arecibo  radar  scatter  diagnostics.  Two  principal 
conclusions  are  offered:  First,  a  wave-particle  scattering  process 
(nonlinear  Landau  damping)  is  responsible  for  the  saturation  of  the 
parametric  instability .  This  theory  provides  a  good  description 
of  the  frequency  spectrum  of  the  scattered  radiation,  especially 
the  plasma  line  portion  (radiation  scattered  with  a  frequency 
shift  almost  equal  to  the  modifier  frequency) .  Second,  a  com¬ 
parison  of  the  intensity  of  the  Arecibo  scatter  diagnostics  with 
theory  suggests  that  the  modifier  power  flux  at  Arecibo  is  below 
threshold  for  the  parametric  instability  but  still  large  enough  to 
produce  enhanced  scattering.  We  recommend  that  ionospheric  modifi¬ 
cation  experiments  be  changed  so  that  the  incident  power  flux  on  the 

2  2 
ionosphere  lies  in  the  100  -  200  ywatt/m  range  instead  of  50  ywatt/m 

2 

(Boulder)  or  20  uwatt/m  (Arecibo)  .  This  increase  in  power  flux  will 
assure  that  short-scale  plasma  turbulence  is  fully  in  the  nonlinear 
range.  We  have  obtained  the  result  for  the  rate  of  absorption  of  the 
modifier  energy  due  to  the  excitation  of  this  turbulence  and  have 
incorporated  this  result  into  the  modifier  propagation  equations  to 
find  the  total  fraction  T  of  modifier  energy  absorbed  due  to  non¬ 
linear  processes.  This  fraction  is  typically  several  tenths. 

One  immediate  extension  of  the  results  should  be  to  the  treat¬ 
ment  of  more  intense  modifier  fields  than  is  currently  possible  in 
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order  to  determine  the  scaling  laws  for  the  turbulence  level  and 
absorption  rate  at  very  high  modifier  fluxes. 

Our  work  on  long- scale  turbulence  concerns  a  theory  for 
artificial  Spread-F,  based  on  a  model  of  a  radiation-induced,  self¬ 
focussing  thermal  instability.  The  basic  idea  is  simple:  A  small 
reduction  in  plasma  density  near  the  reflection  point  for  the 
modifier  radiation  leads  to  focussing  of  the  modifier  radiation 
into  the  region  of  reduced  density.  The  more  intense  radiation 
then  causes  increased  heating  which  via  pressure  gradients  further 
reduces  the  plasma  density.  If  this  model  proves  to  be  correct, 
then  the  cause  of  "artificial  Spread-F"  is  quite  different  from 
natural  Spread-F.  The  predicted  growth  times  are  in  agreement 
with  the  observed  formation  times  (on  the  order  of  1-10  seconds) 
of  magnetic  field  aligned  density  structures.  These  results  can 
profitably  be  extended  to  include  an  analysis  of  shorter  wave¬ 
length  perturbations  than  are  currently  amenable  to  treatment  and 
of  the  nonlinear  behavior  of  the  instability. 

The  controlled  production  of  field-aligned  plasma  density 
striations  allows  one  to  simulate  the  striations  found  in  natural 
ionospheric  disturbances. 
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II.  INTRODUCTION 


The  purpose  of  the  theoretical  effort  at  Princeton  as  sup¬ 
ported  under  ARPA  Contract  No.  F30602-72-C-0053  is  to  obtain  as 
complete  a  theoretical  understanding  as  possible  of  not  only  the 
important  physical  processes  which  lead  to  the  experimentally  ob¬ 
served  effects  in  the  current  ionospheric  modification  program 
(being  undertaken  at  the  Arecibo  Observatory  and  at  Boulder, 
Colorado) ,  but  also  the  scaling  laws  which  predict  the  results  of 
the  future  modification  experiments  (e.g.,  increased  high-frequency 
modifier  power  fluxes) . 

In  the  modification  experiments,  powerful  high-frequency 
radio  transmitters  of  frequency  up  to  the  F  region  cutoff  fre¬ 
quency  fj  are  beamed  at  the  upper  ionosphere.  These  transmis¬ 
sions  are  sufficiently  strong  and  well-collimated  that  even  if 
only  a  moderate  fraction  of  the  incident  energy  were  to  be  ab¬ 
sorbed,  the  energy  balance  of  the  F  region  would  be  signifi¬ 
cantly  changed.  A  review  of  the  experimental  configuration  at 
Boulder  and  many  of  the  results  are  available  in  a  paper  by 
Utlaut  and  Cohen.1 

Extensive  diagnostics  reveal  distinct  collective  phenomena 
on  several  scale  lengths.  The  shortest  scale  phenomena  are  en¬ 
hanced  high-frequency  electron  plasma  oscillations  concomittant 
with  an  increase  in  the  level  of  low-frequency  (on  the  order  of 
the  ion  acoustic  frequency)  oscillations  at  the  same  wavevector. 

At  substantially  longer  wavelengths,  comparable  to  the  ion  Larmor 
radius  r^  (about  5  m)  ,  and  lower  frequencies  (typical  times 


-4- 


for  variation  are  from  1  -  10  seconds),  density  disturbances  with 
wavevectors  nearly  orthogonal  to  the  geomagnetic  field  are  ob¬ 
served. 

2 

In  the  regime  ot  short-scale  turbulence,  our  work  has  shown 
that  the  presence  of  intense  modifier  radiation  will  cause  plasma 
(parametric)  instabilities  in  the  ionosphere.  The  plasma  insta¬ 
bilities  result  in  short-scale  turbulence  involving  plasma  waves 
and  an  increased  absorption  of  the  modifier  power.  "Parametric 
instability"  means  that  electromagnetic  modifier  power  spontane¬ 
ously  converts  itself  into  a  combination  of  electronic  plasma 
waves  and  acoustic  waves.  The  density  fluctuations  associated 
with  the  plasma  waves  and  acoustic  waves  can  scatter  radar  waves 
as  evidenced  by  the  enhanced  signals  seen  by  the  Arecibo  430  MHz 
diagnostics.  Our  theory  has  described  both  the  intensity  and  fre¬ 
quency  spectrum  of  the  Arecibo  radar  scatter  diagnostics.  Two 
principal  conclusions  are  offered:  First,  a  wave-particle  scat¬ 
tering  process  (nonlinear  Landfu  damping)  is  responsible  for  the 
saturation  of  the  parametric  instability.  This  theory  provides 
a  good  description  of  the  frequency  spectrum  of  the  scattered 
radiation,  especially  the  plasma  line  portion  (radiation  scattered 
with  a  frequency  shift  almost  equal  to  the  modifier  frequency) . 
Second,  a  comparison  of  the  intensity  of  the  Arecibo  scatter  di¬ 
agnostics  with  theory  suggests  that  the  modifier  power  flux  at 
Arecibo  is  still  below  threshold  for  the  parametric  instability. 

We  recommend  that  ionospheric  modification  experiments  be  changed 

so  that  the  incident  power  flux  on  the  ionosphere  lies  in  the 

2  2 
100  -  200  uwatt/si  range  instead  of  50  uwatt/n  (Boulder)  ,  or 
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2 

20  ywatt/m  (Arecibo) .  This  increase  in  power  flux  will  assure 
that  short-scale  plasma  turbulence  is  fully  in  the  nonlinear 
range. 

Our  work  on  long-scale  turbulence  concerns  a  theory  for 
artificial  Spread-F,  based  on  a  model  of  a  radiation-induced,  self- 
focusing  thermal  instability.  The  basic  idea  is  simple:  A  small 
reduction  in  plasma  density  near  the  reflection  point  for  the 
modifier  radiation  leads  to  focusing  of  the  modifier  radiation  in¬ 
to  the  region  of  reduced  density.  The  more  intense  radiation  then 
causes  increased  heating  which  via  pressure  gradients  further  re¬ 
duces  the  plasma  density. 

Let  us  now  examine  in  greater  detail  the  status  of  our  under¬ 
standing  of  collective  phenomena  on  the  various  scale  lengths. 
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III.  SHORT-SCALE  TURBULENCE 


The  appearance  of  a  large  enhancement  of  electron  density 
fluctuations  of  large  wavevector  k  ,  comparable  with  an  inverse 
Debye  length  XD  (approximately  1  cm  for  f ^  parameters) ,  i.e., 
,rXD  ~  1  ,  and  of  frequency  close  to  the  modifier  frequency  sug¬ 
gests  that  collective  plasma  wave  phenomena  are  playing  an  impor¬ 
tant  role  in  the  experimental  dynamics. 

At  first  noted  by  Perkins  and  Kaw,  the  energy  density  de¬ 
livered  by  the  Boulder  transmitter  to  the  F  region  at  heights 
near  the  cutoff  exceeds  the  threshold  for  the  excitation  of  a 
parametric  instability  in  which  high-frequency  electron  plasma 
waves  and  low-frequency  ion  acoustic  waves  are  excited  simul¬ 
taneously.  These  instabilities  result  from  the  decay  of  the  modi 
fier  (parametric  pump)  field  into  collective  modes  of  plasma 
oscillation  at  a  rate  faster  than  energy  is  dissipated  from  them 
resulting  in  growth  of  these  modes.  The  wavevector  and  frequency 
sum  rules 


k 


✓wsO 


+  k  . 


/vw  la 


» 


corresponding  to  momentun  and  energy  conservation  respectively, 
must  be  nearly  obeyed  if  the  time  averaged  energy  flow  is  to  be 
nonzero.  Here  the  subscripts  refer  to  the  pump  (o) ,  electron 
plasma  (p) ,  and  ion  acoustic  waves  (ia) ,  respectively.  If  the 
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plus  sign  is  taken  in  the  above  formulae,  yielding  a  plasma  wave 
of  frequency  smaller  in  magnitude  than  the  pump  frequency,  then 
the  direction  of  energy  flow  is  from  the  modifier  into  the  plasma. 
The  opposite  is  tiue  for  the  choice  of  minus  sign.  We  see  then 
that  any  enhanced  level  of  plasma  waves  due  to  this  process  will 
have  a  frequency  less  than  that  of  the  modifier.  Neglecting  dis¬ 
sipation  and  the  geomagnetic  field  for  simplicity,  the  pump  wave 
satisfies  the  dispersion  relation 


.2  2 
+  k  c 
o 


9 


at  each  point  in  space.  The  logitudinal  plasma  waves  satisfy  the 
dispersion  relation 


CO 


2 

P 


[i  +  3<kpy2] 
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where  <o  is  the  electron  plasma  frequency,  and  the  ion  acoustic 
pe 

waves  satisfy  the  relation 


w.  =  ( k A _ )  u)  .  (1 

la  D  pi 


+  3T./Te) 


1/2 


with  w  .  the  ion  plasma  frequency,  and  T  (T . )  the  electron 
pi  e  l 


if  the  sum  rules  are  to  be 


Pi 

(ion)  temperature.  Since  w  >> 

satisfied,  we  must  have  w  =  u>  ,  i.e.,  the  instability  occurs 

op 

where  the  plasma  density  is  near  the  cutoff  density. 

This  instability  has  been  invoked  as  the  cause  of  the  en¬ 
hanced  fluctuation  levels  observed  via  the  backscatter  diagnostics 
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.  3 

at  Boulder.  The  linear  theory  of  this  instability  in  which  one 
derives  the  threshold  (minimum)  level  of  the  pump  field  necessary 
for  growth  of,  and  the  growth  rate  and  frequency  of,  infinitesimal 
perturbations  is  well  known.  In  order  to  apply  a  theory  of  para¬ 
metric  instabilities  to  the  interpretation  of  the  observed  fluctua¬ 
tion  level  the  linear  theory  must  be  supplemented  by  a  nonlinear 
theory  which  takes  account  of  the  finite  amplitude  of  the  excited 
waves  and  the  interactions  between  them  if  the  final  level  of  tur¬ 
bulence  is  to  be  found.  The  formulation  of  equations  describing 
the  nonlinear  evolution  of  the  plasma  wave  turbulence  and  the  so¬ 
lution  of  these  equations  to  determine  the  nonlinear  level  of  the 
spectral  intensity  of  electron  density  fluctuations,  1^  ,  has 
constituted  one  of  the  major  activities  at  Princeton. 

In  order  to  make  progress  in  this  effort,  several  simplifying 
assumptions  were  made  in  the  formulation  of  a  model.  The  principal 
ones  and  their  justification  are  as  follows: 

1.  The  wave  vector  of  the  pump  is  taken  as  zero.  The  di¬ 
pole  approximation  is  a  good  one  because  the  typical  wavevectors 
of  the  electrostatic  plasma  waves  are  larger  in  the  ratio  c/vTe 

3 

(greater  than  10  for  the  F  region)  where  c  is  the  speed  of 
light,  and  vTe  is  the  electron  thermal  velocity.  Thus,  the 
wavevector  of  the  pump  may  be  neglected  in  comparison. 

2.  The  plasma  is  assumed  homogeneous.  Tne  typical  wave- 

1  2 

lengths  of  plasma  turbulence,  on  the  order  of  10  -10  AD  (XD  =  lcm) 
are  small  compared  to  the  density  gradient  scalelength  L  (102  km) 
and  this  makes  the  assumption  plausible.  Indeed,  Perkins  and  Flick^ 
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have  demonstrated  that. the  linear  theory  of  a  homogeneous  plasma 
is  directly  applicable  to  such  a  mildly  inhomogeneous  plasma  as 
the  F  region.  We  return  to  the  question  of  the  effects  of  the 
plasma  inhomogeneity  after  a  discussion  of  the  results  of  the  non¬ 
linear  theory. 

3.  The  pump  amplitude  is  sufficiently  small  so  that  the 

maximum  growth  rate  of  the  unstable  waves,  *  (w  /16) (E2/4imT)  , 

M  pe  o 

does  not  exceed  the  ion  acoustic  frequency,  i.e.,  ym  £  w.  .  Here 
2 

is  the  pump  energy  density,  and  nT  is  the  plasma  particle 
kinetic  energy  density.  For  such  small  pump  intensities,  the  weak 
turbulence  theory'’  of  plasma  has  been  shown  to  be  applicable.^ 

4.  The  electron  and  ion  temperatures  are  comparable 

(Tfi  -  T^) .  This  assumption  (well  satisfied  in  the  ionosphere) 
guarantees  a  large  damping  decrement  for  the  low-frequency  oscil¬ 
lations  (comparable  with  their  oscillation  frequency)  and  makes 
possible  an  explicit  solution  for  these  fluctuations  in  terms  of 
the  high-frequency  plasma  wave  turbulence.  The  problem  then  re¬ 
duces  to  the  simpler  one  of  describing  the  evolution  of  only  the 
plasma  wave  turbulence. 

Under  these  conditions  we  have  shown  that  the  electron  plasma 

wave  spectral  intensity  1^  [the  electric  field  energy  density/ 

3  w 

volume/ (Ak)  ]  satisfies  a  kinetic  equation  of  the  form 


1 

7 


dlk/dt  = 


*k  ^Yk  - 


+  dk'M(k  ,  k')I.  ,]  +  S, 


(1) 
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Here  y^  is  the  linear  growth  rate  due  to  the  presence  of  the  pump, 
v  is  the  damping  decrement  cf  plasma  waves,  M(k  ,  k')  is  the 

c  /w\  /vw 

matrix  element  which  describes  the  nonlinear  interaction,  and 

/vw 

is  the  source  term  which  represents  the  rate  of  spontaneous  emission 
of  plasma  waves. 

Substantial  progress  has  been  made  in  the  solution  of  Eq.  (1) . 
In  general,  1^  is  a  function  of  three  wavevector  components  k.  , 
but  an  examination  of  the  symmetry  properties  of  Eq.  (1)  shows  that 
axisymmetry  of  I.  about  the  direction  of  E  is  preserved  in 

AVvO 

/VW 

time.  We  have  restricted  ourselves  to  a  consideration  of  such  a 

class  of  solutions  making  the  problem  a  two-dimensional  one.  With 

the  use  of  an  approximate  form  of  the  matrix  element  which,  however, 

preserves  the  symmetry  properties  and  magnitude  of  M  ,  we  have 
7 

obtained  an  analytic  solution  to  this  equation  in  two  dimensions. 

We  have  also  analyzed  Eq.  (1)  retaining  the  exact  form  of  M  ,  and 
taking  the  spectral  intensity  as  one-dimensional  with  its  wavevector 
collinear  with  the  pump  field  direction.  In  this  case  we  have 

0 

solved  Eq.  (1)  numerically  and  also  analytically  for  the  time 
asymptotic  solution  in  the  limit  in  which  the  small  source  term 
is  neglected.  The  results  for  the  total  wave  energy  and  nonlinear 
absorption  coefficient  are  in  excellent  agreement  for  the  last  two 
analyses,  and  a  factor  of  three  larger  than  that  obtained  for  the 
analysis  using  an  approximate  form  for  M  We  present  the  results 
of  the  analysis  obtained  by  using  the  exact  form  of  M  below  as 
the  details  of  the  shape  of  the  spectral  intensity  are  more  ac¬ 
curately  obtained  in  this  case.  The  results  are  as  follows: 
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1.  When  the  pump  field  is  switched  on,  the  plasma  wave  energy 
first  increases  exponentially  with  a  growth  rate  equal  to  ym  .  The 
total  wave  energy  reaches  a  greatly  enhanced  level  and  oscillates 
about  a  mean  value  with  an  analytically  predictable  frequency.  These 
oscillations  decay  in  time  as  predicted  by  several  theorems  which 
we  have  proven  regarding  the  general  behavior  of  any  solution  of 
the  kinetic  equation  and  which  demonstrates  that  a  steady  state  can 
be  realized.  This  behavior  is  shown  in  Fig.  1,  where  the  total  plasma 
wave  energy  is  plotted  versus  time  for  a  pump  intensity  about  four 


Fig.  1.  The  evolution  of  the 

total  plasma  wave  energy  density. 

XT  =  I  E^/4irnT  ,  and  ym  is  the 
k 

maximum  value  of  the  linear  growth 
rate.  (722113) 


2.  Figure  2  shows  the  time-averaged  mode  spectrum  for  the 
sample  case  of  Fig.  1.  The  general  characteristics  are  a  distinct 
narrow  spike  near  the  linearly  most  unstable  mode  driven  by  the 
external  pump  field,  followed  by  a  distinctly  separated  spectrum 
of  waves  which  have  been  secondarily  driven  unstable.  In  contrast, 
the  analytic  solutions  obtained  with  an  approximate  form  of  the 
matrix  element  have  shown  a  smooth  spectrum  of  plasma  waves  with 


times  the  threshold  intensity. 


i - '  I  - T - T 
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only  a  single  maximum.  The  primarily  driven  waves  are  separated 
in  frequency  from  the  pump  by  ~a>.  whereas  the  minimum  separation 
from  the  pump  frequency  of  the  secondarily  driven  waves  is  ■“3w^a  . 
The  width  in  frequency  (and  wavevector)  of  the  enhanced  spectral 
intensity  increases  linearly  with  pump  intensity. 


Fig.  2.  The  time-averaged 

spectral  intensity,  1^  ,  of  plasma 

waves  versus  wave  number.  1^  is 

normalized  to  the  plasma  kinetic 

energy  density.  If  the  pump  wore 

a  plasma  wave  of  frequency  tuo  , 

its  wavevector  would  be  k  =  0.17A  1 

o  d 

(722115) 


3.  The  total  wave  energy  increases  as  the  square  of  the  pump 
field  intensity  according  to  the  formula 

IT  =  6.24(Eo/ET)4ve/0.pe  , 

where  the  threshold  electric  field  is 

E^/47rnT  =  16(ve/u>pe)  . 

Here  IT  is  the  total  plasma  wave  electric  field  energy  density 
normalized  to  the  particle  energy  density 


IT  *  I  ( |  E.  |  2/4irnT) 

1  k  K 

Tnis  result  has  been  obtained  both  from  the  numerical  analysis, 
as  is  shown  in  Fig.  1,  and  from  the  analytic  analysis. 

Fig.  3.  The  time-averaged 
energy  density  in  the  self- 
consistent  plasma  waves  versus  the 

a 

pump  field  intensity,  I_  *  £  Ej"/47rnT 

t  .  k 

2  2  K 
and  P  -  e0/et  •  The  solid  line  is 

2 

the  curve  I_  =  6.24P  v  /to  „  . 

T  e  pe 

(722112) 


4.  The  effective  absorption  rate  of  the  pump  v  increases 
linearly  with  the  pump  field  intensity  according  to  the  formula 

v*  =  ve  (Eo/ET)2  *  (2) 

5.  The  result  for  the  local  value  of  the  nonlinear  absorption 
rate  may  be  incorporated  into  the  propagation  equation  for  the 
electromagnetic  wave  in  order  to  determine  the  fraction  of  incident 
energy  absorbed  by  nonlinear  processes.  We  have  done  this  analysis 
using  the  geometric  optics  (WKB)  approximation  for  the  propagation 
of  the  modifier,  assuming  a  horizontally  stratified  ionosphere  and 
a  pump  frequency  less  than  f2.  The  result  (valid  for  I’  <  1)  is 
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r  =  (fq/f*)  m  [  J.  (hoio/c)2/3] 

2 

where  Fq  =  E^c/Qv  =  r.m.s.  power  flux  incident  on  the  ionosphere, 

and  F*  =  nTc(c/Hwo)  .  Here  H  is  the  density  gradient  scale 
2 

length  and  Ev  is  the  vacuum  value  of  the  modifier  intensity. 

Taking  the  following  typical  values  for  the  F  region: 

n  =  4  *105  cm”3  , 

T  -  0.1  eV  , 

H  =  107  cm  , 

7 

a)  =  3*10  /sec  , 
o 

we  obtain  r  »  0.5  for  the  fraction  of  the  incident  energy  absorbed 
by  nonlinear  processes. 

6.  The  question  of  the  angular  dependence  of  the  spectral 
intensity  has  been  resolved  in  part.  The  prediction  of  the  two- 
dimensional  treatment  using  an  approximate  form  of  the  matrix 
element  is  that  the  spectral  intensity  is  very  sharply  peaked  along 
the  direction  of  the  modifier  electric  field.  It  is  expected  that 
a  small  amount  of  inhomogeneity  can  significantly  alter  these  results. 
We  have  done  preliminary  work  along  these  lines  and  have  found  that 
the  extension  of  the  analysis  to  include  a  mildly  inhomogeneous 
plasma  results  in  a  spectral  intensity  with  an  angular  half-width 
on  the  order  of  30°.  Further  work  remains  to  be  done  along  these 
lines. 

There  are  several  areas  in  which  further  research  is  desirable 
at  this  point.  One  is  to  examine  in  more  detail  the  angular  de¬ 
pendence  of  the  spectral  intensity.  Another  is  to  extend  the 
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nonlinear  theory  to  include  more  intense  modifier  fields,  i.e,, 
to  those  which  give  rise  to  linear  growth  rates  somewhat  greater 
than  the  ion  acoustic  frequency.  A  third  area  for  future  re¬ 
search  is  the  inclusion  of  the  magnetic  field  in  the  propagation 
equations  for  the  modifier  and  plasma  waves. 
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IV.  LONG  SCALE  TURBULENCE 

We  have  offered  a  possible  explanation  of  the  formation  of 
magnetic  field  aligned  density  irregularities,  "artificial  Spread 
F,"  as  due  to  self-focussing  of  the  pump  field  into  filaments  of 
varying  intensity  creating  corresponding  density  filaments.  We 
have  formulated  a  model  which  demonstrates  the  physical  processes 
involved  but  which  does  not  treat  the  complex  ray  trajectories 
present  in  the  experiments.  The  principal  results,  those  for  the 
growth  rates  of  such  field-aligned  irregularities,  are  in  reason¬ 
able  agreement  with  the  observed  time  scale  for  density  variation 
on  these  long- length  scales. 

The  mechanism  which  leads  to  instability  may  be  seen  from 
the  following  simple  picture:  In  the  unperturbed  state,  there  is 
no  horizontal  variation  and  the  electromagnetic  flux  reaching  the 
critical  height  is  independent  of  position.  If  we  perturb  this 
state  by  heating  alternate  flux  tubes,  then  the  density  in  these 
tubes  at  the  critical  height  decreases  because  of  plasma  expan¬ 
sion.  The  radiation  is  refracted  into  these  less  dense  tubes 
from  the  adjacent  ones,  causing  further  heating  through  increased 
absorption.  This  is  the  basic  "feedback  mechanism"  which  causes 
the  instability.  The  picture  is  slightly  complicated  by  the  fact 
that,  although  the  density  near  the  critical  height  is  reduced 
causing  ray  convergence,  the  density  further  down  the  flux  tube 
increases  due  to  an  influx  of  particles  from  greater  heights. 

These  two  effects  compete  in  determining  the  net  refraction  of 
the  ray  and  only  if  the  first  dominates  is  there  instability. 
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The  assumptions  which  go  into  the  model  are  as  follows: 

1.  The  ionosphere  is  assumed  horizontally  stratified  and 
the  ambient  magnetic  field  vertical  (in  the  z  direction) . 

2.  The  times  of  interest  are  assumed  sufficiently  small 

and  the  variation  perpendicular  to  the  magnetic  field  sufficiently 
gentle  so  that  the  plasma  motion  across  the  field  may  be  neglected. 

3.  The  modifier  radiation  is  assumed  vertically  incident  and 
totally  absorbed  at  the  critical  density  due  to  the  short-scale 
turbulence  described  above. 

4.  The  energy  in  the  modifier  it  assumed  to  be  preferentially 
delivered  to  high  energy  electrons  (those  in  the  tail  of  the  veloc¬ 
ity  distribution)  which  are  ineffective  in  transferring  their  energy 
to  the  bulk  plasma  until  they  have  traversed  a  mean  free  path.  The 
model  accounts  for  this  by  assuming  absorption  of  all  the  incident 
modifier  energy  flux  at  the  critical  height  and  treating  this  ab¬ 
sorbed  energy  as  a  distributed  source  in  the  electron  energy  equa¬ 
tion.  Since  ion  neutral  collisions  are  felatively  frequent  at  fj 
heights  (vin  -  l/sec) ,  the  ion  and  neutral  temperatures  are  identi¬ 
cal.  Because  the  dominant  energy  loss  mechanism  available  to  the 
electrons,  collision  with  the  ions,  is  slower  lv#i  («,e/iBi)  <<:  vin] 
the  electron  energy  equation  must  be  solved. 

5.  The  electromagnetic  field  is  assumed  to  propagate  in  the 
aeometric  optics  approximation.  This  approximation  limits  our 
treatment  of  perpendicular  wavelengths  to  those  for  which  diffrac¬ 
tion  effects  are  negligible.  This  approximation  is  the  first  to 
break  down  as  the  transverse  wavelength  of  the  irregularity  decreases 
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o.  Since  the  electron  cyclotron  frequency  is  much  less  than 
the  electron  plasma  frequency,  the  magnetic  field  is  neglected  in 
the  ray  equations. 

The  growth  rate  y  was  studied  as  a  functior  of  three  para¬ 
meters:  the  perpendicular  wavevector  of  the  perturbation  p  ,  the 
effective  mean  free  path  of  the  heated  electrons  *  and  the 

sMr;*  of  the  density  profile.  The  latter  variation  was  achieved 
by  choosing  a  parabolic  density  profile  below  the  critical  heigh* 
Zcr  t*1*  form 

n  Z  U-a(Z/Z  cr>I 

for  Z  -  Zcr  .  Here  n£r  is  the  critical  density.  The  values 
a  “  0  and  a  *  °*5  correspond  respectively  to  a  linear  profile 
and  to  absorption  at  maximum  density.  This  profile  is  matched 
along  with  its  first  derivative  at  the  critical  height  to  a  pro¬ 
file  of  the  form 

n  *  (bj  +  b2Z)  exp  (-Z/L) 

for  z  -  Zcr  ,  where  L  ■  ZT/n^g  is  the  barometric  scale  height. 

The  density  profile  is  varied  by  changing  a  ,  keeping  z  fixed 

cr 

at  0.6L.  The  other  important  parameters,  the  electron-ion  and  ion- 
neutral  collision  frequencies  are  assumed  to  vary  with  temperature 
and  density  as  follows: 
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(n/K) 

- TTi  ' 

(T/TQ) i/l 


v 


in 


vir.  eNp  (_2z/L) 


where  K  is  the  peak  density  and  Tq  the  ambient  temperature. 

The  predominant  neutral  is  assumed  to  be  atomic  oxygen.  The  values 
of  =  300/sec  and  v?n  =  4/sec  are  chosen  as  typical . 

The  results  of  tnc  analysis  are  most  conveniently  given  in 
graphical  form  and  are  shown  in  Fig.  4  where  the  growth  rate  is 
plotted  vertically  versus  the  perpendicular  wavevector  of  the  per¬ 
turbation.  The  growth  rate  is  normalised  ''O  the  thermal  diffusion 
2 

time  NL  /R  with  R  =  3 . 16  (NT0/r.'e)  the  electron  thermal  con- 

9 

ductivity  in  the  Z  direction.  The  modifier  energy  density  is 
denoted  by  IQ  and  c  denotes  the  gravitational  acceleration. 

For  reference,  under  typical  conditions  a  value  of  1  for  y  cor¬ 
responds  to  a  growth  rate  of  y  -  i0_2/sec  ,  and  the  dimensionless 
factor  c/v?rg]  lm2/m2] }  in  the  definition  of  p  typically 

equals  10^  . 


Fig.  4.  Plot  of  the  growth 
rate  y  versus  perpendicular 
wavevector  of  the  perturbation 
p,  both  in  dimensionless  units, 
for  several  values  of  the  heated 
electron  mean  free  path  A^  , 
and  the  curvature  a  of  the 
underside  density  profile. 

722121) 
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Problems  for  future  study  include  extension  of  the  analysis 
to  perpendicular  wavelengths  comparable  with  the  ion  Larmor 

radius,  and  formulation  of  estimates  for  final  level  of  turbulence 
due  to  nonlinear  effects. 


t. 
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V.  CONCLUSION 


Our  examination  of  plasma  phenomena  on  both  short  and  long 
length  scales  has  allowed  us  to  obtain  an  understanding  of  many 
aspects  of  the  plasma  physics  of  the  modification  program.  We 
have  developed  a  quantitative  description  of  the  nonlinear  develop¬ 
ment  of  the  short-scale  turbulence  driven  by  the  modifier  near 
the  cutoff  height.  Using  this  description,  we  have  calculated  the 
total  modifier  energy  absorbed  by  nonlinear  processes.  The  ob¬ 
served  field  aligned  density  structures  with  wavelenghts  perpendi¬ 
cular  to  the  magnetic  field  greater  than  an  ion  Larmor  vadius  may 
result  from  a  self-focussing  thermal  instability  induced  by  the 
modifier.  We  have  analyzed  the  linear  theory  of  this  instability 
and  have  obtained  growth  rates  in  agreement  with  the  observed 
formation  time  of  these  irregularities.  Further  work  on  the  non¬ 
linear  evolution  of  this  instability  is  needed.  Additional  work 
is  also  needed  on  the  development  of  a  description  of  plasma 
behavior  in  the  intermediate  scale  length  regime  where  wave lengths 
perpendicular  to  the  magnetic  field  are  comparable  to  an  ion  Larmor 
radius. 
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